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Abstract — We consider the problem of channel estimation 
for amplify-and-forward two-way relays assuming channel reci- 
procity and M-PSK modulation. In an earlier work, a partially- 

] blind maximum-likelihood estimator was derived by treating the 
data as deterministic unknowns. We prove that this estimator 
approaches the true channel with high probability at high signal- 
to-noise ratio (SNR) but is not consistent. We then propose 

] an alternative estimator which is consistent and has similarly 
favorable high SNR performance. We also derive the Cramer- 
Rao bound on the variance of unbiased estimators. 

Index Terms: Amplify-and-forward, Channel estimation, 
Cramer-Rao bound. Two-way relay networks. 



I. Introduction 

Amplify-and-forward (AF) Two-Way Relay Networks 
(TWRNs) m have received a lot of attention for their ability 
to improve the spectral efficiency of bidirectional communi- 
cation. An essential requirement in these networks is reli- 
able channel knowledge at both ends of the communication 
link. The channel estimation problem for AF two-way relays 
has been considered in several recent works Q-IH. These 
works developed training-based algorithms which estimate 
the channel using pilot symbols known to both terminals. 
Unfortunately, the required training overhead may significantly 
reduce the spectral efficiency of the system. It is therefore 
worthwhile to develop reliable partially or completely blind 
channel estimation approaches that minimize this loss of 
spectral efficiency. 

In |5l and ||6l, the authors developed partially -blind 
maximum-likelihood (ML)-based channel estimation algo- 
rithms for TWRNs assuming that the channels are flat fading 
and that the two terminals employ A/-PSK modulation. In 
both works, the data symbols were treated as deterministic 
unknowns and pilots were only needed to resolve the phase 
ambiguity. The difference between f5'| and Q is that the 
channels were assumed to be reciprocal in Q but nonrecip- 
rocal in j6Jl. It was proved in IS) that the ML estimator for 
nonreciprocal channels is consistent and that it approaches the 
true channel with high probability at high SNR for M > 2. 

In this work, we focus on reciprocal channels. We analyze 
the asymptotic high SNR behavior of the ML estimator 
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developed in ||5] and investigate its consistency. We show that 
for AI > 2, the ML estimator approaches the true channel 
with high probability at high SNR. However, we also prove 
that it is not consistent. No such analysis was provided in ||5l. 

As an alternative to the ML estimator, we propose to 
estimate the channel by minimizing the sample variance of 
the envelope of the received signal after self-interference 
cancellation. This criterion was used by the ML estimator 
developed for nonreciprocal channels in [6:|. In this work, we 
investigate its application to reciprocal channels. We term the 
proposed estimator as the minimum sample envelope variance 
(MSEV) estimator. The asymptotic behavior of the estimator 
established in ||6l still holds under channel reciprocity. Thus, 
the MSEV estimator is consistent and approaches the true 
channel with high probability at high SNR for M > 2. 

As a reference, we also derive two Cramer-Rao bounds 
(CRBs) on the variance of unbiased estimators under chan- 
nel reciprocity. Monte-Carlo simulations are then used to 
obtain the mean-squared error (MSE) of the two estimators, 
demonstrating that the proposed estimator outperforms the ML 
estimator In summary, the main contributions of this work are: 
(0 analysis of the high SNR performance of the ML estimator 
for reciprocal channels; (//) investigation of the consistency of 
the ML estimator for reciprocal channels; (;//) application of 
the MSEV criterion to reciprocal channels; (iv) derivation of 
two CRBs on the variance of unbiased estimators for reciprocal 
channels. 

The remainder of this correspondence is organized as fol- 
lows. In Section Ull we present our system model. In Sec- 
tion Unl we present the ML estimator of Q and the proposed 
estimator In Section IIVI we analyze the high SNR behavior 
and the consistency of the two estimators. The CRBs are 
derived in Section [V] Our simulation results are shown in 
Section [Vl] Finally, our conclusions are in Section IVIII 
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Fig. 1 . The two-way relay network with two source terminals and one relay, 
assuming channel reciprocity. 
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II. System Model 

We consider the half-duplex TWRN with two source termi- 
nals, 7i and T2, and a single relay, TZ, shown in Fig. [1] The 
network operates in quasi-static flat-fading channel conditions. 
Each data transmission period is divided into two phases |[T| . 
In the first phase, 7i and T2 simultaneously transmit to TZ 
the M-PSK data symbols h = ^/l\e^'*'^ and ^2 = V7^e•''^^ 
respectively, where Pi and P2 are the transmission powers, and 
(f>i and (j)2 are chosen from the set S^j = {{2£— 1)tt/M,£ ~ 
1, . . . , M}. The received signal at the relay during this phase 
is given by r = hti + gt2 + n where h and g are the complex 
coefficients representing the flat-fading channels Ti TZ 
and T2 ^ TZ, respectively, and n is the zero-mean complex 
additive white Gaussian noise (AWGN) with variance cr^. In 
the second phase, the relay broadcasts the amplified signal Ar, 
wher43 A = \J p^^^^^ ^- Without loss of generality, we will 
consider estimation at 71. Assuming channel reciprocity, the 
received signal at 7i during the second phase is 



N 



\h\ = ^//p- J2 - Aatii\, and that a is given by 



i=i 



z = Aati + Abt2 + Ahr 



(1) 



where a = h?, h = gh, and i] is the zero-mean complex 
AWGN term at Ti with variance cr^. Since, A/-PSK modula- 
tion is assumed, it is sufficient for detection to know a and 



III. Channel Estimation Algorithms 



Let the vectors z = [zi...zn]'^, ti = [in, . . . , iijv]"^, 
*2 = [^21, • ■ • ,^2^]^, n = [ni,...nN], and r] = [771, . . . , Ty^v] 
denote the received vector at 7i, the transmitted symbol 
vectors of 7i and T2, the noise vector at TZ and the noise 
vector at 7i during N successive transmissions, respectively. 
We note that the sample size N is constrained by the coherence 
time of the channel during which the channel parameters a and 
b remain constant. 



A. Maximum Likelihood Estimation 

We begin by briefly presenting the ML estimator for 
reciprocal channels derived in |l5]. To avoid dealing with 
a complicated likelihood function, the transmitted symbols 
t2i,i = 1,. . ■ ,N are treated as deterministic unknowns. We 
also ignore the finite alphabet constraint that restricts the 
phases 4>2i ~ ^t2i, i = l...,iV to the set Sm. The ML 
approach can be used to blindly estimate the parameter a and 
the sums -0^ = 02i + i'bi i = I, ■ ■ . , N. However, it cannot be 
used to obtain separate estimates of 0b and (j)2i, i — 1, . . . ,N. 
Hence, a small number of pilots will be needed to isolate 
an estimate of (pb, and because of this we will refer to the 
ML estimator as partially blind. The unknown parameters are 
collected in the vector = [a, . . . , V^jv]"^. Let a, \b\, 

ipi be the ML estimates of a, \b\, and tpi, i = 1,...,A^, 
respectively. It is shown in IS) that tpi = Z{zi — Aatu), 

'This choice of the amplification factor maintains an average power of Pr 
at the relay over the long term (2)- 



N 



a = argmm ■ 



N 



■J2 \ Autiil - if J2 - Autik 

i=l \ fe=l 

(j2(A2|u| + l) 

+ TV log (^2 + 1) |. 

(2) 

We show in Section |IV] that the estimator in (|2]l is not 
consistent. We next propose an alternative criterion which 
results in a consistent estimator of a. 

B. Proposed MSEV Algorithm 

Let Zi[u) = {zi — Autxi) , i ~ 1, . . . ,N he the "cleaned" 
versions of the received signal samples after self-interference 
has been removed using the complex value u as an estimate 
of a. The signals Zi{u), i = 1,...,N sue independently 
generated realizations of the random variable z{u) ^ A{a — 
u)ti + Abt2 + Ahn + ?/. The quantit}0 
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WNiu) 



JV / TV \- 

1=1 \ k=l / 

represents the sample variance of the envelope of z{u). One 
would intuitively expect the actual variance of \z{u)\ to be 
smallest in the absence of self-interference (i.e., at u = a). We 
therefore propose the MSEV estimator of a which is given by 

1 ^ / 1 ^ \' 

i=l \ k=l ) 

The MSEV estimator has the same form as the ML estimator 
of the self-interference channel for the nonreciprocal case 
in [51. However, it does not retain the same maximum- 
Ukelihood interpretation when channel reciprocity is assumed. 

The solutions for (O and ^ may be obtained using iterative 
methods such as the steepest-descent algorithm or the quasi- 
Newton algorithm Q. Since the objective functions in (|2| 
and (|4|i are nonconvex, the performance of such methods 
will depend on the availability of good initial estimates. A 
simple way to obtain an initial estimate is the sample average 
estimator used in [6] and given by = jj^j^ X^i^i ^h-^j- 
In our experimental studies, we will use the steepest descent 
algorithm to solve ^ and (|4|. As we shall see, the resulting 
MSE for both estimators is almost identical to that obtained 
when the solutions are acquired using grid search. In the 
next section, we compare the asymptotic behavior of the two 
estimators. 

IV. Asymptotic Behavior Analysis 

Our analysis focuses on the estimation of a since |&| is not 
required for detection. We study the asymptotic behavior of 
the two estimators by considering their consistency and their 
behavior at high transmit SNR. The transmit SNR is defined 

^Wn{u) is a scaled version of the numerator of the first term in {2). 
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7 — As we shall see, both d and d^, approach a as 
7 ^ oo with high probability for AI > 3. Moreover, the 
MSEV estimator is consistent while the ML estimator is not. 
We start our analysis with the MSEV estimator 

A. High SNR Behavior of the MSEV Estimator 

It can be easily verified that the high SNR behavior of 
the ML estimator under nonreciprocal channels in |[6l (see 
Lemma 1 and Lemma 2 therein) is not affected when channel 
reciprocity is assumed. Hence, letting G^{u) ^ ^WV(u), 
and assuming that Pi ~ aP2 and = f3{Pi + P2) for some 
a, /? > 0, the following theorem holds. 

Theorem 1 For fixed N, and assuming M-PSK modula- 
tion, lim G^{u) = +00 for all u ^ a with probability 

2 ^^^'^ [M - 1). For u^a, lim G^{a) < +00. 



1 - 



7—)- 00 



The above theorem guarantees that for M > 3 and for 
sufficiently large N, the MSEV estimator will approach the 
true channel with high probability as SNR increases. 

B. High SNR Behavior of the ML Estimator 

Denoting by A-y(u) the objective function of the ML esti- 
mator in (|2]i, we have that 

A^{u) = NG^iu)/iA^\u\ + 1) + A^log {A^\u\ + l) 
iV((l + a)7 + l)G'7(u) 



I3{1 + a)-f\u\ + {1 + a)j + 1 
'I3{l + ah\u\ 



(5) 



A^log 



Hence, 



(1 + a)7+ 1 
N{l + a) 



lim A^(u) =—, s. . , : — 

7->oo ^(1 + Q,)|y| _)_ (1 _)_ Q,) 7-^00 

+ N\og{l3\u\ + l). 



lim G~f{u) 



(6) 



Thus, by Theorem [T] lim A™ (a) < +00 and lim K^{u) = 

7— )-oo 7— foo 

+00 for all a with probability 1 - ( ^g) ^ (M- 1), i.e., 
the ML estimator exhibits similar high SNR behavior to the 
MSEV estimator. 

C. Consistency of the MSEV Estimator 

The proof for the consistency of the ML estimator under 
nonreciprocal channels in Q is not affected when channel 
reciprocity is assumed. Since Wn{u) is the sample variance 
of it converges in probability to the variance of 

which is given by IS): 

W(m) =A^\a- MpPi + A2|5|2P2 + a^{AM + 1) 

^a\AM + ^) f . f S (7) 



. k=0 



where Li/2{x) is the Laguerre polynomial lH with parameter 
1/2, Ir{-) is the Modified Bessel Function of the First Kind 
of order t, and 

1 



cT^A^\a\^ + l) 



{A'\v\^Pi+A'\b\'P2 



(8) 



2A^\v\\b\y/PiP2Cos{Zv ~ (j)b + 2iTklMy 



The variance 'W{u) has a unique global minimum at u = a, 
and the estimator is consistent when the channel parameters 
a, b belong to compact sets. 

D. Consistency of the ML Estimator 

Let Y/v(u) be the objective function of the ML estimator 
in (|2]i scaled by ^y^J. Thus, 

Wn{u) N 



Yn{u) 



-log {A^ 



(t2(A2|u| + 1) N- 
As iV — > 00, Y/v(u) converges in probability to 

W{u) 



y{u) 



0-2(A2|u| + l) 



log(A2 



(9) 



(10) 



To find out whether y{u) has an extremum at w = a, we an- 
alyze the behavior of the partial derivative^ J^sr [u) = 

and 3>q( u) — ■ Because of the symmetry of y{u) with 

respect to 5R{u} and it is sufficient to consider y^{u). 

We have 

A^^i{u} 



u\{A^\u\ + \) 

1 



1 - 



W{u) 
a^{A^\u\ 
dW{u) 



1) 



(11) 



ct2(A2|u| + 1) d^{u}' 



Moreover, it can be shown that 3^^/"? 



0. Therefore, 



For the factor 



A^^{a} 



W{a) 



(12) 



1 - 



(proved in Appendix 



W(a) 



a2(A2|a| + 1) 
we have the following lemma 



Lemma 1 For any a G 



1 _ W(a) 
^ a^A^\a\ + l) 



> 0. 



Hence, yu{a) is zero only when 5R{a} = 0. Similarly, y<i{a) 
is zero only when 3{a} is zero. Since y{u) is differentiable 
at u = a (for a 7^ 0) and u = a is not a boundary point, this 
implies that y{u) does not have an extremum at u = a for 
a 7^ 0. Hence, the ML estimator is not consistent |9|. 

The inconsistency of the ML estimator should not come as 
a surprise since the data symbols are treated as deterministic 
unknowns. Due to this assumption, the received samples 
are not identically distributed and the number of estimated 
parameters is not fixed but grows linearly with the number 
of samples. Hence, the well-known sufficient conditions for 
the consistency of ML estimators (see ifTOl Theorem 10.1.6]) 
are not satisfied. The unknown data symbols are referred to 
in the literature as incidental parameters, while the channel 
parameters which are the same for all samples are called 
structural parameters ifTTI . lfT2l . It is well known that the 
asymptotic properties of ML estimators, such as consistency, 
do not necessarily hold in the presence of incidental parame- 
ters IfTTI . In fact, the ML estimator of structural parameters can 
be inconsistent even when a consistent estimator exists 1 12]. 



^In this work, we use } and Sfj } to denote the real and imaginary 
parts of a complex number, respectively. 
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V. Cramer-Rao Bounds 

We derive two Cramer-Rao bounds for the estimation 
problem under consideration. The first bound is derived 
by treating the phases Vi i ■ • ■ : i^N as deterministic un- 
knowns. The vector of unknown real parameters is 0/; = 
[3?{a}, 3{a}, V-i, • ■ ■ , fpNV and the PDF of z is 

f{z \en)= ^ 



exp 



7r(T2(A2|a| 



~ Aatii 



(13) 



(72(A2|a| + l) 



Let I{0r) be the corresponding Fisher information matrix 
(FIM). We have 



91og/(^ I eH)\fdlogf{z I Or 



B 

C 



(14) 



where 



A = 



A^n{a}'3{a}N 



B' = 



2A2 



ct2 (^2|a| _,_!) 



A'^SR{a}Q{a}Af 
|a|^(A^|a| + l)^ 



3{6*tljvt2w} ^{b*tiNt*2j^} 



(16) 



and 

C = 



a2 (A2|a| + i; 



[2A2iVP2, 2A2|&|2p2, . . . , 2A2|6|2p2] 

(17) 



Assuming I{9r) is invertible, the CRB for the estimation 
of a is given by the sum of the first two diagonal entries 



in the inverse of 1(0;?), i.e0, CRBa 



CRBa = tr ( [A- BC^B' 



[I (0_r)]22- Let D be the 2 X 2 top left submatrix of 
I^^{9r). Using the Schur-complement property, we have that 

D= (A-BC^^B'^) , i.e.. 



(18) 



Because the symbols ti are known and the symbols t2 are 
treated as deterministic unknowns, CRBa is a function of ti 
and t2, and it thus applies for the particular realizations of ti 
and t2 under consideration. 

Another variant of the CRB commonly used in the pres- 
ence of random nuisance parameters is the Modified CRB 
(MCRB) in). In deriving the MCRB, we compute a mod- 
ified FIM (MFIM) by first extracting the submatrix of the 
conventional FIM which corresponds only to the nonrandom 
parameters (Jftja}, 5{a} and \b\) and then obtaining the expec- 
tation of this submatrix with respect to the random nuisance 
parameters (f2 in our case). Let 6' = [^{a},^{a}, \b\]'^, and 



A. 



*The notation [A]ij is used to refer to tlie (i,j)th element of tiie matrix 



let I{9') be the corresponding 3x3 submatrix of I{6 r). Then 
the MFIM is given by (19). The MCRB on a is the sum of 
the first two diagonal entries in the inverse of I,n{d'), i.e., 

MCRBa^[i;„\e')U + [i-,\9%2 

_ 4a^Pi{A^\a\ + 1)2 + g4^2(^2|^| ^ (20) 
iNA^Pf{A'^\a\ + l) + 2Na^A^Pi ' 

While the above bound has the advantage of a simple closed 
form, it is not as tight as CRBa. 

VI. Simulation Results 

In this section, we compare the MSB performance of the 
ML estimator and the MSEV estimator using Monte-Carlo 
simulations. Our results are generated assuming AI = 4 
(QPSK) and Pr = Pi = P2, and are averaged over the 
same set of 100 realizations of the channel parameters g, h 
which are independently generated from the complex Gaussian 
distribution with mean zero and variance 1. The minimizers 
of the objective functions in (|2]i and (|4|i are obtained using the 
steepest descent algorithm. The initial points are chosen using 
the sample average estimator, while the step size is chosen 
using backtracking line search |7:|. As a reference, we also 
show the MSB performance for the two estimators when the 
solutions for Q and ^ are obtained using grid search with 
a step size of lO^'^. We also show the bounds CRBa and 
MCRBa, where CRBa is averaged over many realizations 
of ti and t2- 

Fig. I2] shows the average MSB of the two estimators vs. 
SNR for N = 100. The MSBV estimator outperforms the ML 
estimator and the performance gap is most significant at low 
to medium SNR. At high SNR, both estimators approach the 
bound CRBa- We can also see from Fig. |2] that, for both 
algorithms, the MSB performance that results when steepest 
descent is employed is almost identical to that when grid 
search is used. 

The bar plots in Fig. [3] show for both estimators the average 
number of steepest descent iterations required for convergence 
and the average number of line search iterations to find the 
step size for a single steepest descent iteration, respectively. As 
we can see from Figs. |2l [3] the steepest descent algorithm is 
a reliable and low complexity method for solving (|2]i and (|4|. 
Moreover, the MSBV estimator is more efficient than the ML 
estimator since on average it requires fewer steepest descent 
iterations to converge and fewer line search iterations to obtain 
the step size. 

Fig. |4] shows the average MSB of the two estimators vs. N 
for an SNR of 15dB. The MSBV estimator has a superior MSB 
performance which improves steadily as N increases. The gap 
between the MSB performances of the two estimators becomes 
more significant as the sample size increases, in accordance 
with the fact that the MSBV estimator is consistent while the 
ML estimator is not. 

VII. Conclusions 

In this work, we compared two partially-blind channel 
estimation algorithms for two-way relay networks assuming 
channel reciprocity and A/-PSK modulation. The first was 
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7,„(0') " EW)} = 



<J-^(A^\a\+l) ^ |a|^(A^|a| + l)^ 
A''5R{a}a{a}Af 
|aP(A^|a| + l)^ 





y4*5R{a}Q{a}Ar 
|a|^(A^|a| + l)^ 








2A'^NP2 
T^(A^|a|+l) 



(19) 



the ML estimator obtained by treating the data symbols as 
deterministic unknowns H, and the second was the MSEV 
estimator which minimized the sample variance of the enve- 
lope of the received signal after self-interference cancellation. 
We showed that both estimators approach the true channel 
with high probability at asymptotically high SNR and that 
the MSEV estimator is consistent while the ML estimator is 
not. We also derived two CRBs on the variance of unbiased 
estimators. Monte-Carlo simulations were used to compare 
the MSE performance of the two estimators, showing that 
the MSEV estimator performs better than the ML estimator 
and that the steepest descent algorithm can be used to provide 
accurate low-complexity implementations for both estimators. 



Appendix A 



We have 

W(a) 



1 - 



(T2(A2|a| + l) 



1/2 



A^\h\^P2 



A^\b\^P2 

o-2(A2|a| + 1) 



(21) 



Let Q{x) = J (-^1/2 i^^)) ^ 2;, it is sufficient to show that 
Q{x) > for a; > 0. In [9, Appendix E], it is shown that 
Q{x) is strictly decreasing for a; > 0. Using this fact, we 
will establish that Q{x) > for a; > by showing that 
lim Q{x) = -i. We can expand Q{x) as 



Q{x) = -e-^ 



{l + xyio{- +2xil + x)Io{- h{- 



(22) 



Moreover, for large arguments, we have the following asymp- 
totic expansion [Sl 



Ir{x) 



(£-!)(£ -9) 



, (23) 



/2ttx L 8a; 2!(8x)2 

where e = 4t2. Using the above expansion, we obtain the 
following approximations for large x 

.2 X 5 

4 + 8' 
x 1 
2 + 4' 

.2 X 3 

4 ~ 8' 

Substituting these approximations into ( l22l i. we have 



^«-2.(l + x)7„(f)A(f 



(24) 
(25) 
(26) 
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lim Q{x) = i, which completes the proof. 



Fig. 2. Average MSE of the ML estimator and the MSEV estimator and the 
bounds CRBa and MCRBa vs transmit SNR for N = 100, and A/ = 4 
(QPSK). For each estimator, two MSE plots are shown: one for the grid 
search (GS)-based chaimel estimate and one for the steepest descent (SD)- 
based estimate. 
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Fig. 3. Average number of steepest descent iterations (a), and average number 
of line search iterations for a single steepest descent iteration (b) required by 
the ML and MSEV estimators for different SNR values (N = 100, M = 4). 




Fig. 4. Average MSE of the ML estimator and the MSEV estimator, and 
the bounds CRBa and MCRB a vs sample size lor an SNR of 15 dB, and 
M = 4 (QPSK). For each estimator, two MSE plots are shown: one for 
the grid search (GS)-based channel estimate and one for the steepest descent 
(SD)-based estimate. 



